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1. VISHNU hybrid model: 

Viscous hydrodynamics is a useful tool to describe the expansion of the QGP fireball created in relativistic heavy 
ion collisions. Previous viscous hydrodynamic simulations reveal that elliptic flow vi and triangular flow vz 
are very sensitive to the shear viscosity. Even the "minimal" specific shear viscosity given by the KSS bound 
r\js = 1/4-7T [1] can lead to a significant suppression of «2 [2, 3] and U3 [4, 5]. In principle, this allows for an 
extraction of the QGP viscosity from viscous hydrodynamic calculations by tuning rj/s to fit the experimental 
i>2 and U3 data at RHIC and LHC. In practice, this requires excellent control of other inputs in the theoretical 
calculations, which include the initial conditions, the equation of state (EOS), hadronic dissipative effects, and a 
proper description of the chemical and thermal freeze-out in the hadronic stage [6]. 

For a more realistic description of the late hadronic stage, we developed the hybrid code VISHNU [7], which couples 
viscous hydrodynamics for the QGP expansion with a hadron cascade model for the late hadronic evolution at 
a switching temperature T sw . VISHNU consists of three parts: VISH2+1 solves the Viscous Israel-Stewart Hy- 
drodynamic equations in 2+1 space-time dimensions, with an assumption of longitudinal boost invariance [2], 
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H20 [7] converts the hydrodynamic output into particles profiles for further hadron cascade (UrQMD) propagation 
through a monte-carlo event generator based on the modified Cooper-Frye formula, and UrQMD solves the Boltz- 
mann equation for a variety of hadron species with flavor-dependent cross-sections [8]. By describing the hadronic 
rescattering and freeze-out procedure microscopically, VISHNU improves earlier, purely hydrodynamic models and 
eliminates the additional adjustable parameters required for the transport and freeze-out characteristics of the 
hadron phase, making it possible for a reliable extraction of the QGP viscosity from the experimental data. 
The default switching temperature in VISHNU is set to T sw — T c h = 165 MeV, which is close to the QGP phase 
transition temperature and the chemical freeze-out temperature at RHIC energies. Using a partial chemical 
equilibrium EOS as input, we investigated whether the microscopic Boltzmann approach can be emulated by 
the macroscopic hydrodynamic approach by varying the switching temperature T 3W in VISHNU simulation, and 
found a strong T^-dependence of the elliptic flow (with a constant rj/s as input). We then extracted a tempera- 
ture dependent effective shear viscosity (r]/s) e ^ (T) for UrQMD and found that pure viscous hydrodynamics with 
(v/ s ) e (T) as input could nicely fit the p-r-spectra and differential elliptic flow V2(pr) from VISHNU. However, 
the extracted (T) strongly depends on the pre-hydrodynamic history, particulary the chosen value of the 

QGP shear viscosity. It therefore does not represent the intrinsic transport properties of the hadronic matter, 
but a parameter that reflects some memories of the QGP transport properties. We conclude that there exists no 
switching window below T c h in VISHNU, where viscous hydrodynamics can replace the hadron cascade and remain 
predictive power; Both components of the hybrid model are needed for a quantitative and realistic description of 
the dynamical evolution of the QGP fireball created in relativistic heavy ion collisions [7]. 

2. The QGP shear viscosity at RHIC and LHC energies 

During the evolution of the QGP fireball, the shear viscosity controls the efficiency to convert the initial spacial ec- 
centricity e x to the fluid momentum anisotropies e p . However, the distribution of the final momentum anisotropics 
to the flow patterns «2(pt) for different hadron species strongly depends on the chemical composition and the 
accumulated radial flow in the late hadronic stage. In Ref. [9, 10], we proposed to use the integrated elliptic flow 
t>2 h for all charged hadrons to extract the QGP shear viscosity since it is directly related to the fluid momentum 
anisotropies e p and less sensitive to other details in VISHNU calculation: such as the chosen bulk viscosity, the 
detailed form of the viscous Sf corrections, etc. Ref. [9] also shows that the theoretical V2 h /e x curves as a function 
of multiplicity density per overlap area dN c h/(dyS) are universal, which depend only on the specific QGP shear 
viscosity (tj/s)qgp, but not on the details of the initial conditions. Furthermore, pre-equilibrium flow and bulk 
viscosity only slightly affect these theoretical curves, which are at the order of 5% or below 5%. It thus preferable 
to extract the QGP viscosity from a comparison between the theoretical and experimental V2 h /e x — dN c n/(dyS) 
curves, as shown in Fig. 1. The colored lines with symbols are the VISHNU results with different {t]/s)qgp as 
input. Left and right panels correspond to two different smoothed initial conditions, obtained from averaging a 
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large number of fluctuating initial entropy densities (given by MC-Glauber or MC-KLN models) by aligning the 
participant plane for each event. For the experimental curves, V2 and dN c h/dy are experiment measurements, 
while e and S are theoretical inputs obtained from the MC-Glauber and MC-KLN models. This leads to the 
differences in magnitude and slight changes in slope for the two experimental curves shown in Fig.l left and right. 
As a result, the extracted value of (j)/s)qgp from these two panels changes by a factor of 2 mainly due to the 
different e x from MC-KLN and MC-Glauber. Recent event-by-event viscous hydrodynamic simulations showed 
that the triangular flow V3 is also sensitive to the QGP shear viscosity [4, 5] which suggests that, in the near 
future, the combined analysis of V2 and V3 from VISHNU could yield an even more precise value of (ti/s)qgp and 
may reduce the ambiguities in initial conditions from the hydrodynamic side. At this moment, we take the current 
uncertainties from the initial conditions and conclude from Fig.l that 1 < 4ir(r)/ s)qgp < 2.5 for the temperature 
region T c < T < 2T C probed at RHIC. 

With (rj/s)QGP — (l/47r) for MC-Glauber and {t]/s)qgp — (2/47r) for MC-KLN extracted from thepr integrated 
V2 for all charged hadrons, VISHNU yields an excellent description of the p-r-spectra and differential elliptic flow 
V2(pr) (for pt < 2GeV) for all charged hadrons and identified hadrons at different centrality bins measured in 
200 A GeV Au+Au collisions at RHIC [10], showing the power of the newly developed hybrid model and the 
robustness of the above extraction of the QGP shear viscosity. 

The successful fit of the RHIC low pr data allows for constrained LHC predictions from VISHNU. After 
extrapolating our calculation to the Pb+Pb collisions at LHC and comparing with recent flow data from the 
ALICE Collaboration, we found that the (t]/s)qgp extracted from RHIC slightly over-predicts the LHC data. 
After increasing {ti/s)qgp from 2/(47r) to 2.5/(47r) (for MC-KLN), VISHNU gives a better description of the 
integrated and differential elliptic flow data for all charged hadrons at LHC [11]. A further simulation from 
VISHNU also yields a reasonable description of the V2(pr) for identified hadrons such as pion, kaon and protons 
at different centralities [12]. These results suggest that the averaged specific QGP shear viscosity (or the 
assumed constant (v/ s )qgp) at LHC energies is slightly larger than at RHIC energies. However, a further study 
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assuming different temperature dependencies for (t)/s)\qgp shows that one cannot uniquely constrain the form of 
(v/ s )\qgp(T) by fitting the spectra and «a alone. Based on our current understanding, the question on whether 
the QGP fluid is more viscous or more perfect in the temperature regime reached by LHC energies is still open [11]. 
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